Diffusive Lorentz gases and multibaker maps 
are compatible with irreversible thermodynamics 
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We show that simple diffusive systems, such as the Lorentz gas and multibaker maps are perfectly 
compatible with the laws of irreversible thermodynamics, despite the fact that the moving particles, 
or their equivalents, in these models do not interact with each other, and that the dynamics takes 
place in low-dimensional phase spaces. The interaction of moving particles with scatterers provides 
the dynamical mechanism responsible for an approach to equilibrium, under appropriate conditions. 
This analysis provides a refutation of the criticisms expressed recently by Cohen and Rondoni 
[Physica A 306 (2002) 117-128]. 
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I. INTRODUCTION 



PACS: 05.20.-y; 05.45.+b; 05.60.+w; 05.70.Ln 



For nearly one hundred years, the Lorentz gas has been investigated as a model of diffusive transport of light tracer 
particles among heavier particles [1-5]. On the collisional time scale of the fast tracer particles, the heavy particles 
can be treated immobile, which simplifies the dynamics. Due to this simplification, the Lorentz gas plays a privileged 
^ and important role in the development of transport and kinetic theories. In the last decades, several versions of the 
Lorentz gas model have been mathematically studied in detail and the existence of a well-defined diffusion coefficient 
has been proved rigorously under certain conditions [6,7]. Moreover, the Einstein relation between the coefficients of 
-H ' electrical conductivity and diffusion has also been proved for a Gaussian thermostated Lorentz gas in the presence of 
J> | an external field [8,9]. 

Ten years ago, one of us introduced the multibaker map as a simple model of deterministic diffusion [10]. This 
map is constructed by simplifying the Birkhoff-Poincare map of the hard-disk Lorentz gas. The multibaker map is a 
spatial extension of the baker map described by Arnold and Avez [11] and previously by Hopf [12] who also described 
a spatially extended baker-type model identical, in many respects, to the multibaker map. 1 The advantage of the 
, multibaker map for a study of deterministic diffusion is that, like the baker map, it is exactly solvable and can be 
used to test specific hypotheses in the context of kinetic and transport theories. Several versions of the multibaker 
| map have been proposed and investigated as models of different transport processes [10,14-27]. 
. — ■ Recently Cohen and Rondoni have written a series of papers [28-30] where, among other things, the physical 
relevance of the Lorentz gas and multibaker map in irreversible thermodynamics is challenged. Their main criticism 
• • is that the Lorentz gas and multibaker map "represent noninteracting particle systems" so that they "do not possess 
the crucial property of local thermodynamic equilibrium" . Hence "these models are not suited for a derivation of 
irreversible thermodynamics" (Ref. [29] p. 117) where, in their view, "the local physical density gradient dn/dx acts 
like a real force which does work on the particles themselves". In contrast, as claimed, "in the Lorentz slab there is 
no dissipation, no work done on the moving particles, hence no entropy production" (Ref. [29] p. 127). 

Our purpose in this paper is to refute these criticisms. We show that the views of Cohen and Rondoni concerning 
both the status of the density gradient as a real force performing work on the particles and the physical unsoundness 
of the Lorentz gas and multibaker map from the point of view of irreversible thermodynamics, are wrong. Matyas, 
Tel, and Vollmer have recently responded to criticisms of Cohen and Rondoni, directed at their work as well as ours, 
and expressed views along lines similar to those given here [31]. 

The plan of the paper is the following. In Section II, we give a brief review of the irreversible thermodynamics of 
diffusion in binary mixtures and clarify the status of the thermodynamic forces present. In the following Sections III 
and IV, we show that diffusion in the Lorentz gas and in the multibaker map is well described by, and compatible 
with irreversible thermodynamics, and that a local equilibrium indeed establishes itself in these systems. 



1 The baker map was introduced in the thirties in the work by Seidel [13] who seems to have been inspired by Birkhoff. 



1 



II. THE IRREVERSIBLE THERMODYNAMICS OF MATTER EXCHANGES 



In this section, we provide a first response to the claim that thermodynamic forces do work on the particles 
themselves and that this is a sine qua non condition for irreversible entropy production. We briefly review the 
irreversible thermodynamics of isothermal diffusion in multi-component mixtures and show how it applies to Lorentz 
gases, considered as binary mixtures of light and infinitely massive particles. Moreover, we show that the entropy 
production of diffusion does not vanish in binary mixtures of particles with an extreme mass ratio. 



A. The Second Law in the presence of matter exchanges 

Classical irreversible thermodynamics is a local theory. It rests on the sole assumption that the densities of entropy 
s, internal energy e, and particle numbers n k are locally linked through a relation having the same structure as the 
relation linking the global analogs of these quantities in equilibrium: 

s(f,t)=s[e{f,t),{n k {f,t)}} . (1) 

The differentiable function s has the same properties as the entropy 

S = S(E,V,{N k }) = Vs(y,{^ , (2) 

of a system at the internal energy E and containing N k particles of substance k in a volume V. Here, the identity of the 
two differentiable forms allows one to interpret the derivatives of s with respect to e and n k as being, respectively, the 
local reciprocal temperature and minus the chemical potential of substance k multiplied by the reciprocal temperature. 

Starting from Eq. (1) one can derive a balance equation of entropy by differentiating both sides with respect to 
time and by using the balance equations of mass, momentum, and energy of macroscopic physics. One obtains then 
in the balance equation a source term - the entropy production - in the form 

CTirr = J2 J {a) X^ > , (3) 

a 

where is the thermodynamic flux associated with process a and is the conjugate thermodynamic force. In 
a given system all these processes need not coexist. When they do they are, typically, coupled as long as certain 
symmetry requirements are satisfied. 

The status of Eq. (1) from the standpoint of microscopic theory has been the subject of many investigations. 
Within the context of kinetic theory of gases, Prigogine [32] showed that a sufficient condition of validity is that the 
velocity distribution must remain close to the Maxwellian. This guarantees, at the same time, the validity of linear 
phcnomcnological relations linking the fluxes to the forces, a property that is less general than Eq. (3). However, 
Eq. (1) may still be correct under circumstances where the local equilibrium distribution is not a near-Maxwcllian 
distribution in the energy [33]. The near-Maxwellization is sufficient for the validity of (1) but docs not constitute 
a necessary condition. As a matter of fact, the main physics behind Eq. (1) is that after a short initial relaxation 
period the space-time variation of the probability distribution descriptive of the nonequilibrium system at hand can 
be expressed entirely in terms of the variation of the macroscopic (hydrodynamic) fields. 

If a system exchanges matter dAf k as well as heat dQ with its surroundings, the Second Law reads [34] 



AS = S Il -S I > h^ + > S k dM k , (4) 



where AS is the total change of entropy in the system, not including the surroundings, S k is the partial molar entropy 
of substance k, and Af k denotes the number of moles of this substance. In a closed system that does not exchange 
matter with its surrounding (i.e., dAf k = 0) but only heat, we recover the well-known inequality 



AS = S u -Si> 



which implies that if no heat is transferred to a closed system, there may still be a transient irreversible increase of 
entropy. On the other hand, in an open system, matter exchange contributes to entropy production beside processes of 
heat exchange. In this regard, we mention that the study of processes of exchange of matter requires the consideration 
of chemiostats (particle reservoirs) instead of thermostats (heat reservoirs). 
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B. Irreversible thermodynamics of diffusive processes in mixtures 



Transport by diffusion is a process in which there is an exchange of matter from one macroscopic part of a system 
to another. This irreversible process can take place in a system in which no other irreversible process occurs. We then 
have a process where there is an exchange of matter contributing to entropy production without further contribution, 
for instance, due to temperature or macroscopic velocity gradients. Processes of exchange of matter are very common 
in chemical thermodynamics [35]. In particular, the process of mixing of two distinct chemical species is well known to 
contribute to entropy production. In such processes, spatial inhomogeneities of chemical concentrations exist across 
the system before the thermodynamic equilibrium is reached. The diffusion of tracer particles in a fluid of other 
particles is an example of such an irreversible process with exchange of matter. 

A fluid composed of tracer particles moving among other particles can be considered as a binary mixture. It is well 
known that an isothermal c-component mixture at rest (i.e., with a zero macroscopic velocity field) can only sustain 
(c — 1) independent diffusive processes [34,36]. In a binary mixture, we should thus expect a single diffusive process. 
It is a process of mutual diffusion between the tracer particles and the other particles. The position of the tracer 
particles is defined with respect to the center of mass of the total system. 

The thermodynamics of a c-componcnt mixture is based on the Gibbs relation 

c 

dE = T dS - P dV + VkdN kl (6) 

k=l 

where T is the temperature, P the pressure, and \ik the chemical potential of particles k. If we define the densities of 
energy, entropy and particles k respectively as [cf. Eq. (2)] 

E S N k 

e = y , s = y , and n k = — , (7) 



we obtain the local Gibbs relation [37] 



de = T ds + ^ jj,k dnk ■ (8) 



fc=i 



The extensivity of the energy, entropy, and particle numbers together with the local Gibbs relation (8) imply the 
Gibbs-Duhcm relation 

c 

s dT - dP + Y n k dfik = , (9) 
fe=i 

which shows that, for an isothermal and isobaric system (dT = dP = 0), the variations of the chemical potentials are 
interconnected by 

c 

^n fc V^ fe = 0. (10) 

fe=i 

On the other hand, in absence of chemical reactions, local conservation of particles holds, yielding 

d t n k + V • (n k u + Jkj = , (11) 

where Jk is the current density of particles k with respect to the local barycentric velocity u of the fluid, which is the 
velocity of the center of mass of a fluid element. The local conservation of the total mass implies that the currents 
are related by 

c 

Y, m kJk = 0, (12) 
fc=i 

where nik denotes the mass of one particle of component k. 

As a consequence of the constraints (10) and (12) only (c — 1) gradients and currents are independent because 



3 



c-1 



fe=i nc 

In an isothermal c-component mixture at rest (u = 0), the entropy production is therefore given by [36] 

c c-1 

fc=i fc=i 

with the thermodynamic forces 

X fe = -vf-g^vf. (16) 
;=i 

In a binary mixture (c = 2), the entropy production contains only the contribution from the single diffusive process 
the system can sustain 2 

(c=2) ^ = -(i + ^U.?" ( 17 ) 

V m 2 n 2 J T 

If the tracer particles of mass mi are much lighter than the other particles of mass 777,2 and, moreover, if the tracer 
particles are dilute 

777l Tl\ < 777 2 77 2 , (18) 

the entropy production reduces to the simple expression 

CT irr = - Ji • , (19) 

corresponding to the simple diffusion of the tracer particles. The heavy particles fix the center of mass with respect 
to which the transport is defined and they play the role of scatterers for the lighter particles. The thermodynamic 
flux of this diffusive process can thus be identified as the current of the tracer particles 

J 1D) = Ji , (20) 

while the associated thermodynamic force is given in terms of the gradient of the chemical potential of the tracer 
particles as [38] 

X<P) = -V^ . (21) 

Moreover, for dilute tracer particles, the chemical potential depends logarithmically on the tracer density n\ according 
to 

f = *Bln^, (22) 

where n\ is a reference density [35]. Hence, in a dilute system where the temperature, T, is constant, the thermody- 
namic force becomes 



2 Equation (17) is given by de Groot and Mazur [36] as Eq. (94) of Chapter XI in the form <Ti rr = — -^Ji ■ V^r where 0,2 is 
the mass fraction of substance 2 which is indeed given by 0,2 = m<2n<2l{min\ + 7712772). In the Lorentz-gas limit (7771/7772) — > 0, 
the mass fraction of substance 2 equals 0,2 = 1. We notice that we are here using quantities per unit volume while de Groot 
and Mazur use quantities per unit mass in Ref. [36]. Both derivations lead to the same entropy production a per unit volume, 
as it should. 
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(23) 



which is independent of the temperature. The thermodynamic force is simply given in terms of the gradient of tracer 
density, and drives the irreversible process of matter exchange at the macroscopic level. This irreversible process occurs 
even if the tracer particles are so dilute that they never meet each other. The point is that the thermodynamic force 
of diffusion is a concept associated with a collection of particles, in contrast to the microscopic forces which directly 
act on the particles themselves. As we show in Sec. Ill A, such a force is indeed manifested as a Newtonian force at 
the level of momentum balance of the light particles but the type of balance one refers to here pertains to the average 
momentum taken over the probability distribution of the particles, as usual in fluid mechanics and in macroscopic 
physics in general. It may be worth recalling that in some molecular dynamics-inspired mechanisms of thermostating 
designed to sustain nonequilibrium steady states with shear or heat flow, the noncquilibrium constraint (an analog of 
the local thermodynamic force considered here) does modify the dynamics of the individual particles [39,40]. In our 
view, this ingenious scheme, which stimulated many developments, is to be looked at as a prototypical model along 
with other thermostating mechanisms [41-46]. In this connection we mention that nonequilibrium steady states with 
processes of matter exchange would require nonequilibrium constraints on particle numbers, that is, chemiostating 
mechanisms. 

In the linear range of irreversible processes, the flux is proportional to the thermodynamic force. We shall write this 
proportionality in a form displaying Fick's diffusion coefficient V rather than Onsager phenomenological coefficient, 

J 1 ^) = -V Vni , (24) 

since V varies weakly in a wide range of values of n\. Then the contribution of diffusion to the irreversible entropy 
production is given by 

a irr = J<°) • X<P) =k B V > . (25) 

m 

The diffusion coefficient T> can be positive and finite even in the case of dilute tracer particles in a fluid. Since each 
tracer particle undergoes a random walk due to collisions with the surrounding fluid particles, a statistical collection 
of them appear to be driven from the regions where they are more concentrated toward the regions where their 
concentration is lower. The statistical character of the random walk is at the origin of the irreversible production 
of entropy. Thus, the absence of mutual interaction between the tracer particles does not prevent the existence of a 
thermodynamic force contributing to entropy production. 

In conclusion, a careful application of the notion of thermodynamic force in the case of diffusion, as explained 
above, shows that the Cohcn-Rondoni argument that, in the case of diffusion, the thermodynamic force should act 
on the particles themselves [29] is incorrect. It results from a failure to appreciate the fact that the thermodynamic 
force of diffusion has a statistical origin and should not be conceived as a microscopic force which acts on the 
particles themselves. Indeed, the tracer particles can be so dilute that they never come within interaction range 
among themselves although the thermodynamic force of diffusion is still non-vanishing and contributing to entropy 
production. Moreover, the derivation here above shows that the entropy production does not vanish in the Lorentz-gas 
limit (mi/m 2 ) — ► where the heavy particles are immobile. 

In the next section, we show starting from a microscopic description that diffusion in the Lorentz gas is indeed well 
described by irreversible thermodynamics, and that the entropy produced corresponds to physical work that might 
be done on a moving piston, for example. 

III. THE IRREVERSIBLE THERMODYNAMICS OF THE LORENTZ GAS 

In this section, we show, contrary to the claims by Cohen and Rondoni [29], that the Lorentz gas can approach 
an equilibrium state through a local equilibrium and that a thermodynamic force as well as a friction force do exist 
in the Lorentz gas. Furthermore, we show that diffusion in the Lorentz gas yields a positive entropy production and 
that such an entropy production has a clear physical implication. 

A. ii-theorem for the random Lorentz gas 

The Lorentz gas model was introduced by Lorentz in 1905 as a model for transport of electrons in a solid [1]. It 
consists of a collection of fixed scatterers, and a collection of moving particles which do not interact with each other, 
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but only with the fixed scatterers according to a well-defined potential energy function. Today, the Lorentz gas is 
often used as a model of diffusion of light particles among heavy ones [2-5]. The heavy particles are immobile on 
the time scale of motion of the light particles, and the light particles are considered as not interacting with each 
other. At each elastic collision between light and heavy particles, energy is nearly conserved although momentum is 
not. Therefore, the randomization of the velocity direction is very fast albeit the randomization of energy occurs on 
a much longer time scale. This decoupling of time scales does not preclude the existence of a well-defined diffusion 
coefficient on each energy shell. 

To set the stage we first analyze the transport properties and the thermodynamics of the Lorentz gas from the 
standpoint of kinetic theory. To this end we consider the random Lorentz gas with dilute scatterers, a version of the 
Lorentz model in which the diffusion coefficient has been evaluated some time ago using kinetic theoretic methods 
[2,5]. More specifically we consider a two-dimensional dilute random Lorentz gas. The scatterers are assumed to 
be hard disks of radius a which are uniformly distributed on the plane with the density n^. The motion of tracer 
particles is described by a linear Boltzmann equation, also known as the Boltzmann-Lorentz equation [2,47] 



dtf + v-Vf 



n 2 av 



+7T 



dip' 



sin ■ 



\f(?,<p')-m<p)] 



(26) 



for the number density f(r, Lp) of tracer particles at position r with velocity angle <p, where v = v(cos<p, siny>) is the 
velocity of the particles, v their speed which is conserved as well as their kinetic energy e = \rri\v 2 . 
The Boltzmann-Lorentz equation (26) obeys a iJ-theorem based on the local entropy density 



sir) 



+ 7T g , 

dip f(f, Lp) In — - 



f(r,<p) ' 



(27) 



where f° is a constant, and e = 2.718 . . . denotes the Naperian base. Indeed, the balance equation for this entropy 
density is given by 



d t s 



-V • J s + a u 



(28) 



with the entropy current 



/ + 7T 
d<pv f(r,ip) In-— - 
f{r, 



(29) 



and the local entropy production 



CTirr = T &B ™2 



avj 



dLp dLp' 



sin ■ 



Lp -LD 1 



/(*W) 

f{r,<P) 



> . 



(30) 



We emphasize that this 77-theorem holds although all particles have the same energy e — \ra\v 1 . The mechanism 
at the origin of this if-theorem is the randomization of the velocity angle, which establishes a local equilibrium in 
velocity direction. We notice that the existence of such a local equilibrium has already been pointed out by Lebowitz 
and Spohn in their work on heat conduction in the random Lorentz gas [48]. 

One can expand the density function as a Fourier series in the velocity angle Lp as 



(31) 



As a consequence of the Boltzmann-Lorentz equation (26) the Fourier components fk obey the following coupled 
differential equations 



at, d x -id y d x + \d y 
dt fk + v — y - fk-i + v f k+1 = 



8fc 2 



n 2 av f k ■ 



(32) 



2 ' ~ " ' " 2 ■'^ ± I -4k 2 
For particle distributions close to the microcanonical equilibrium, we keep only the first Fourier components so that 



2 

ni(r) + ^ « ' -Mf) 



(33) 
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where the density of tracer particles and their current are given by 



/+TT 
/(r, if) dip , (34) 

Ji(r) = f + vf(r,<p)dip. (35) 

J — 7T 

Substituting expansion (33) in the Boltzmann-Lorentz equation (26) and retaining the dominant terms, we get the 
balance equations for the number of tracer particles and for their current as 

d t m = -V ■ Ji , (36) 
8 v 2 - 

d t Ji ~ -- n 2 av Ji - — Vni . (37) 

Over a time scale longer than the relaxation time t = 8n 3 2av , the current is driven by the gradient of density and we 
obtain from Eq. (37) Fick's law on the energy shell as 

Ji ~ -V Vni . (38) 

After substituting in Eq. (36) we obtain the diffusion equation 

dtn\{f) = V V 2 ni(r) , (39) 

with the energy-dependent diffusion coefficient 

16 ani 



where v = \J1tjm\ is the speed of particles. Under these conditions, the entropy production (30) is straightforwardly 
shown to reduce to 

a irr ~ k B JL ~ k B V > , (41) 

Vni n\ 

which is the expression expected from irreversible thermodynamics, here shown to hold on a single energy shell. 

We next show that the thermodynamic force of diffusion (23) appears explicitly in the balance equation for the 
momentum of the gas of tracer particles. We notice that the momentum density g\ of the tracer particles is directly 
proportional to the particle current J\ according to 

<?i =m 1 J 1 . (42) 

The balance equation (37) can thus be interpreted as the balance equation for the momentum of the gas of tracer 
particles. Moreover, we can introduce the partial mass density of the tracer particles as p\ — mini and their drift 
velocity as 

ui = — = — . (43) 
Pi m 

This drift velocity can be considered as the velocity of an clement of the gas of tracer particles. Since Eq. (37) gives 
the balance equation for momentum, we can now derive the time evolution of an element of the tracer gas as 

dui _ _ ^± , . . . 

Pi —jj- - ~P Pi ui -e Vni (44) 

-ffriction -^thcrmo. 

where we have neglected terms which are cubic in the gradients of the density. We notice that, to the order we 
keep, neglecting the cubic terms allows us to replace the partial time derivative in (37) by a total time derivative in 
(44). The constant (3 = ^n 2 av is the friction coefficient, while e is the kinetic energy of the particles. Eq. (44) can 
be interpreted as Newton's equation for the element of the tracer gas. The left-hand side contains the mass of the 
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gas element multiplied by its acceleration while the different Newtonian forces (per unit volume) acting on this gas 
element can be identified in the two terms of the right-hand side. 

The first term is the Newtonian force of friction. The presence of this term contradicts the claim by Cohen and 
Rondoni that "fluxes of particles are moving frictionless" in the Lorentz gas (Ref. [29], p. 121). On the contrary, the 
infinitely heavy particles which are the scatterers produce a friction on the tracer gas because of the interaction of 
each tracer particle with the scatterers. We notice that conservation of total energy is here reduced to the property 
that the microscopic kinetic energy of the individual particles, e = em, is locally conserved according to 

d t e + V-(eu 1 )=0. (45) 

Friction terms will be manifested when e is decomposed into a macroscopic kinetic energy part, n\v\/2, and an 
internal energy part. This is exactly along the same lines as the derivation of the balance equations of macroscopic 
physics. 

The second term in the right-hand side of Eq. (44) involves the gradient of the density and can thus be identified as 
the Newtonian force driven by the thermodynamic force of diffusion (23). Indeed, we have a proportionality between 
both given by 

F thcrmo . = 211 JO) . (46) 
k B 

This result shows that, in the Lorentz gas, the thermodynamic force acts as a "real force" albeit on the statistical 
collection of particles which compose the gas element, even if the particles have no mutual interaction and, furthermore, 
conserve their energy during their whole motion. This is very similar to the role of pressure gradient in the momentum 
balance of a volume element of a fluid, which is present even in the limit of noninteracting particles. It is clear that 
this thermodynamic force does not act on the individual particles themselves. The concept of thermodynamic force 
applies to the statistical level of description but not to the microscopic one. 

In the stationary state, the two terms in the right-hand side of Eq. (44) cancel, and one has 

ui ^ Vni . (47) 

The entropy production (41) of diffusion can then be expressed in terms of the dissipation arising from the friction 

<7 irr ~ fc B ^ ul . (48) 

In Ref. [29], Cohen and Rondoni claim that there is no thermodynamic force for diffusion in the Lorentz gas, because 
there is no local equilibrium state described by a Maxwell-Boltzmann distribution with space and time dependent 
macroscopic quantities - the local density and temperature, and the local mean velocity. Here, we have demonstrated 
that despite its simplicity, the random Lorentz gas has all the ingredients needed for a description of its transport 
processes by irreversible thermodynamics. Thus the claims to the contrary by these authors appear to be mistaken. 
In our opinion they have not made a sufficiently careful distinction between the microscopic and macroscopic levels 
of description of the system, and they have not appreciated the fact that the approach to equilibrium in a Lorentz 
gas with a small density gradient also proceeds through a state of local equilibrium. Here the local equilibrium is not 
a state of near-Maxwellian velocity distribution but, rather, a state which is nearly uniform in the velocity direction, 
with a local density, even though all of the moving particles have the same energy. This local equilibrium is at the 
very origin of diffusion itself and, thus, of irreversible entropy production on each energy shell. 

Finally, we notice that the random Lorentz gas is a useful model for the electronic conductivity at low temperature 
in a solid with impurities, and leads to a determination of the so-called residual resistivity. In this case, it is known 
that the conductivity is related to the diffusion coefficient and that both should be evaluated on the energy shell 
corresponding to the Fermi energy of the electron gas [49] . Diffusion as well as entropy production on a single energy 
shell have therefore experimental relevance. 

B. The periodic Lorentz gases as models of diffusion 

We now proceed to the extension of the previous refutation of the criticisms by Cohen and Rondoni to the case of 
the periodic Lorentz gases where diffusion was rigorously proved. Using the methods of dynamical systems theory, 
which are free of the statistical hypotheses underlying kinetic theory, we show in the present subsection that, for these 
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Lorentz gases, local equilibrium in the form discussed in Sec. Ill A and Fick's law also hold on each energy shell and, 
in the following subsection, that this results in a positive entropy production. 

In 1980, Bunimovich and Sinai proved that the two-dimensional hard-disk periodic Lorentz gas has a positive and 
finite diffusion coefficient if the horizon is finite, i.e., if no trajectory exists running through the lattice of hard disks 
without elastic collision [6] . A similar result has been proved by Knauf for the two-dimensional periodic Lorentz gas 
composed of scatterers with Yukawa potentials under the condition that the energy of the moving particle is large 
enough [7]. In addition, these Lorentz gases have been proved to be ergodic and mixing. The mixing property holds 
on each energy shell, and leads to the establishment of a local equilibrium in velocity direction. As in the case of the 
random Lorentz gas discussed above, this local equilibrium may be described by a spatial and velocity distribution 
which is nearly uniform in the directions of the velocity of the moving particles, and the spatial density varies only 
over distances large compared to the mean free path of the particles. This local equilibrium is sufficient for a transport 
by diffusion to occur, compatible with the laws of irreversible thermodynamics. 

Mathematical proofs of the existence of diffusion in periodic Lorentz gases given by Bunimovich and Sinai and 
Knauf [6,7], as well as extensions to higher orders in density gradients by Chernov and Dettmann [50], all imply the 
existence of a well behaved diffusion equation for these models. That is, on large spatial scales, compared to the 
mean free path of a particle, the diffusion occurs on each energy shell so that the density f(r, e) of light particles with 
energy e at position r obeys the diffusion equation 

d t f(r,e)~V(e) V 2 f(f,e) , (49) 

where T>(e) is the diffusion coefficient on the shell of energy e, as proved in Ref. [6]. 

The diffusion equation (49) can be rewritten in the form of a local conservation law for the density f(r, e) as 

8 t f + V • J f = , (50) 

with the current 

J f ^-V(e)Vf(f,e). (51) 

This result, which was directly derived from the microscopic dynamics in Ref. [51], shows that Fick's law holds on each 
energy shell in the Lorentz gas. This is a consequence of the property of local equilibrium in the velocity direction 
on each energy shell, and responds to the doubts expressed by Cohen and Rondoni [29] about the existence of a local 
equilibrium state and a corresponding Fick's law in the periodic Lorentz gases. 

At each elastic collision the energy e of the particle is conserved so that the particle keeps its energy. Particles 
with different energies can thus be distinguished. Hence, the gas can be envisaged as a mixture of different gases 
with particles having different energies. These different gases ignore each other and each of them is the stage of an 
irreversible process of diffusion according to Eq. (49). In such a gas, the number of light particles per unit volume 
around the spatial point f is defined as 

r-OG 

m(r)= / f(f,e)de. (52) 
Jo 

We can consider separately the gases of particles with different energies. If we remove from the system all the particles 
except those having their energy between e and e + Ae (with Ae arbitrarily small), the particle density is given by 

m(r)~/(f,e) Ae. (53) 

As a consequence of the diffusion equation (49), the density (53) obeys the diffusion equation 

d t m(r) =V VV(r) , (54) 

with the diffusion coefficient V = T>(e). We notice that Eq. (54) as well as Eq. (49) are also found in the random 
Lorentz gas with the diffusion coefficient (40) according to the kinetic theory of the previous subsection. 

The diffusion equation (54) describes processes in which the initial state is created by a spatial disturbance in 
the particle density n\. Thereafter, the system relaxes back to a spatially uniform state according to Eq. (54). 
In this spatially uniform state, all the particles still have the same energy e while their velocity angle is uniformly 
distributed in the interval [0, 27r[ and their position is uniformly distributed in the whole space accessible to the 
particles. This uniform state is a microcanonical state for each one of the particles. This state plays the role of the 
state of thermodynamic equilibrium in the Lorentz gas. 
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C. The thermodynamics of the Lorentz gas 



Here, we show by using the methods of thermodynamics that diffusion in the Lorentz gases has a positive entropy 
production on each energy shell. We first establish the equilibrium thermodynamics of the Lorentz gas in the micro- 
canonical ensemble and subsequently derive the balance equation for the entropy density in order to obtain the local 
entropy production. 

The equilibrium thermodynamics of the hard-ball Lorentz gas is established as follows. Let us consider a gas of 
N± light particles without mutual interaction and moving in a spatial domain Q bounded internally by the hard balls 
and externally by a wall enclosing the Lorentz gas in a square, say. The volume of the domain Q is V = Vol(Q). We 
notice that this volume is related to the physical volume V by 

V = V - N 2 v 2 , (55) 

where N 2 is the number of hard balls and v 2 is the volume of one hard ball. The phase space of this gas is 

M = |(fi,pi, ...jfjVi.PivJIrj e Q,pj 6 R d ,xi < x 2 < ■ ■ ■ < x Nl } , (56) 

where fj = (xj,yj, Zj) and the last conditions provide the ordering of indices required in the case of identical particles. 
The phase-space volume element is 



dX = dn 



dr Nl dpi ■ ■■ dp Nl . 



(57) 



The system being ergodic, the equilibrium invariant distribution is given by a microcanonical distribution for a system 
of Ni independent particles of kinetic energy between e and e + Ae (with Ae arbitrarily small) : 



F(X) = 



I 



r 



(27rTOie) 2 Ae 
otherwise , 



for e < 



2mi 



< e + Ae 



(58) 



which is normalized according to 



f F(X) dX = 1 



The entropy can be calculated as Gibbs' coarse-grained entropy [52] 

S = -k B ^2pi\npi , 



(59) 



(60) 



by coarse graining the phase space into cells of size AX = A Nld r A Nld p — (ArAp) Nld . This is the microscopic 
definition of the entropy which has been adopted in Refs. [54,55] for nonequilibrium distribution. The probability for 
the system to belong to the i th cell is 



F dX ~ Fi AX , 



th coll 



so that 



~ k B ( 



In- 



\ F{X)AX I 

A straightforward calculation of the entropy (62) for the equilibrium distribution (58) leads to the formula 

Ve (27rmie)^ Ae 



S ~ k B Ni In 



JVi L (f ) e {ArAp) d ' 



(61) 



(62) 



(63) 



in the limit N\,V — » oo, keeping the ratio Ni/V constant (where, again, c = 2.718 . . . denotes the Naperian base). 
The energy of the gas is 
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E = N 1 e 



(64) 



We notice that S is a well-defined differentiable function. As a consequence, differentiating both the formula (63) for 
entropy and the energy equation of state (64) with respect to the variables V, Aq, N 2 , and e, we obtain the following 
identity which is known as the Gibbs relation [36] 



with the inverse temperature 



the equation of state 



d S = ^dE+^dV-^dN 1 -^dN 2 , (65) 



T = — ~ ' (66) 



Nl fcB (67) 



r v - n 2 v 2 ' 



and the chemical potentials given by 



T (V -N 2 v 2 ){2ir mi e)i Ae 



for the light particles and by 

^ 2 Nik^n, 



r v - n 2 v 2 



(69) 



for the hard balls. It follows from Eqs. (68) and (69) that ii\/T and H2/T, which are related to the thermodynamic 
force of diffusion, are well defined. Notice that the inverse temperature happens to vanish in the two-dimensional case 
[cf. Eq. (66)], while the ratios of pressure and chemical potential to temperature are non-zero. 
If we introduce the intensive quantities defined by Eq. (7), we obtain the other Gibbs relation 



where the entropy density is given by 



with the reference density 



ci- 1 T(f) e(ArAp) d 
In Eq. (70), the energy density is given by 



ds = ^ de — dni — ^ dn 2 , (70) 



s = k B ni In ° 2 Hl , (71) 



o _ (1 - n 2 v 2 ) (27rTOie) 2 Ae 



the light-particle chemical potential by 



the hard-ball chemical potential by 



and the inverse temperature by Eq. (66). 



e = ni e , (73) 



Y = fcBln i' (74) 



[i 2 nifc B ^2 

T=T^r 2 > (75) 



11 



The relation (70) holds for the system at equilibrium. Its local equilibrium analog would amount to a supposition 
that the density ri\ of light particles presents small variations over large spatial scales while the energy e and the 
hard-ball density n 2 remain constant, and the entropy density s varies in such a way that the Gibbs relation (70) 
remains satisfied locally. The relaxation to such a state of local equilibrium finds its full justification in the complete 
ab initio derivation of entropy production given in Ref. [55]. The condition of validity is that the randomization of 
the velocity direction is faster than the relaxation of spatial inhomogencitics of the particle densities on each energy 
shell. The randomization time is given by the intercollisional time t co ii. ~ i/v which is estimated by the ratio of the 
mean free path I to the speed v — y/2e/m\. The relaxation time is given by the diffusion time idiff. ~ L 2 /T> where L 
is the wavelength characterizing the spatial inhomogeneities. Since T> ~ v£ , the condition idiff. S> t co h. thus requires 
that the wavelength of the particle density is larger than the mean free path: L>£ This condition is independent 
of the speed so that it holds uniformly on all the energy shells. 

For a noncquilibrium situation with a local variation of the density n\ of light particles, we can now obtain the 
balance equation for the entropy density (71) on the basis of the knowledge that the particle density n\ obeys the 
diffusion equation (54) in the Lorentz gas. The density n 2 of hard balls is constant in time and uniform in space. 
The entropy balance equation is obtained by taking the time derivative of the entropy density (71), replacing dtn\ 
by the right-hand side of Eq. (54), and splitting the result into a divergence term and a rest term which can then be 
identified as the local source of entropy, i.e. the irreversible entropy production. Performing this calculation yields 
the same balance equation as (28) but here with the entropy current 

J s = -Vk B I In 02 - Ul j Vni = —V Vs , (76) 



and the local irreversible entropy production 



a iTI = k B V > , (77) 

ni 

which is always nonnegative by the Second Law of thermodynamics. The conclusion is here that the standard 
thermodynamic reasoning to derive entropy production [36] can be developed in the Lorentz gas and yields the same 
result as Eq. (41) obtained by kinetic theory in the case of the random Lorentz gas. The present derivation is more 
general in the sense that it also applies to the periodic Lorentz gas. 

We now consider a gas of tracer particles moving in the same system as above, though with a nontrivial distribution 
of energy f(f, e) obeying the diffusion equation (49). As mentioned above, such a gas can be considered as a mixture 
of different gases with particles having their energy defined within intervals Ae. The local entropy density of this gas 
is given by the sum of the local entropies (71)-(72) for the particles with different energies. Since the density of the 
particles of energy e is given by ni(r) ~ /(r, e)Ae, we get the local entropy density of the whole gas as 

f°° , , £ ( —* \ 1 (1 - n 2 v 2 )e(27rTOie)i 
sif) = ale k B f(r, e) In — . (78) 

Jo f{r,e) T (f) e{ArAp) d 

By performing a similar calculation as above but using the diffusion equation (49) , we can obtain a balance equation 
similar to Eq. (28) now for the full entropy (78). In this case, the local irreversible entropy production in the Lorentz 
gas is given by 

r - 1 2 

roc V/(r, e) 

a irr (r) = / dek B V(e) 1 — ^ > . (79) 
Jo f{r,e) 

which is one of the central results of this paper. This result shows that the Lorentz gas with a distribution of energy 
can be considered as a mixture of different gases distinguished by the energy of their particles and that diffusion on 
each energy shell of the gas separately contributes to a positive entropy production. The main point is that, already in 
the microcanonical ensemble, the entropy production is positive with the typical form (77) expected from irreversible 
thermodynamics. The origin of this entropy production holds in the randomization of the velocity angle due to mixing 
in the Lorentz gas, which leads to the local equilibrium in the velocity direction. 

In Refs. [17,18,53-55], it was shown that a similar calculation can be carried out ab initio at the microscopic level 
of description, which leads to the very same expression (77) that we have here derived using the local Gibbs relation 
(70) and the diffusion equation (49). 
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In conclusion, equations (41), (77) or (79) clearly demonstrate that entropy production does not vanish in the 
Lorentz gas, contrary to the claim by Cohen and Rondoni in Ref. [29] p. 127 that there is no entropy production in 
the Lorentz gas. 

We notice that our result (79) applies not only to the Lorentz gas in which the mass ratio m 2 /mi is infinite but 
also to a binary mixture where the mass ratio m 2 /mi is large but finite. In this latter case, there still exists a 
mechanism which thermalizes and drives the velocity distribution toward the Maxwell-Boltzmann distribution even if 
this mechanism is slow. Under this assumption, the local distribution of energy can be supposed to be factorized into 
a local particle density m (r) multiplied by a uniform equilibrium Boltzmann energy distribution of given temperature 
T as 

-(e) e -e/k B T 

(80) 

where 77(e) is the density of states at energy e. Consequently, the local entropy density (78) takes the form 

, ( —\ 1 (! - n 2 v 2 )ci +1 (2nmik B T)i 

s(r) = Kb nAr) In - - - — — — -, . (81) 

w w m(f) (ArAp) d v ; 

expected from the Sackur- Tetrode formula where Ar Ap = 2-kTi [56] . This result shows that our definition of entropy 
corresponds to the well-known physical entropy. Moreover, the entropy production (79) reduces to the form 

cr irr ~ k B V ^ Ul) > , (82) 

m 

with the mean diffusion coefficient 

_ J °°,(e)e-'/^P( e ) de 

J™v(e)e- £/kBT de " { ' 

With the assumption (80), Eqs. (82) and (83) thus give the values of entropy production and diffusion coefficient 
for the diffusion of light tracer particles in a dilute binary mixture at the temperature T, as can be derived from 
the results by Chapman and Cowling in Ref. [4]. We can therefore conclude that our whole calculation is consistent 
with common knowledge on irreversible thermodynamics and the kinetic theory of binary mixtures of light and heavy 
particles. 



D. Thermodynamic equivalence with a process producing work 

We have already shown with Eqs. (44) and (46) that, in the random Lorentz gas, the thermodynamic force appears 
as a Newtonian force acting on statistical collections of particles. In the present subsection, our aim is to give a 
general argument demonstrating that the entropy production of diffusion in the Lorentz gases has a physical meaning, 
since it can be associated with the performance of work. We respond in this way to the doubts expressed by Cohen 
and Rondoni about a possible "confusion between information theoretic and real entropy" (Ref. [29] p. 118), and 
the allusion that entropy production in the Lorentz gas could only refer to an abstract information theoretic concept 
without physical consequence (Ref. [29] p. 118 & 128). 

Let us consider a slab of the Lorentz gas of width I between two large cubic reservoirs of the same volume V = L d 
containing respectively and N-r light particles so that the particle densities in the reservoirs are respectively 
«l = Nl/V 7 and 7Jr = Nr/V (The subscript ! L' refers to the left-hand side and 'R' to the right-hand side, see 
Fig. 1). We here assume that the velocity distribution of the particles is initially taken to be a Maxwell-Boltzmann 
distribution at the same temperature T across the whole system, i.e., the distribution function has the form (80) at 
the initial time with n\{r) taking the constant values and 7Jr in the reservoirs and a linear interpolating profile 
inbetween. Since the diffusion process conserves the energy of each particle, the temperature, T, at the beginning is 
equal to the temperature after equilibrium is reached and there is no heat produced. However, there is a current of 
particles J*- -* from the high-density reservoir to the low-density one because of the difference of concentrations across 
the Lorentz slab. The current of particles decreases the difference of densities between both reservoirs. As long as 
there is a difference of densities, the system is out of equilibrium, and thermodynamic equilibrium is reached when 
Ul(£ = 00) = nn(t = 00) = n oq = (TVl + Nr)/(2V), and the current vanishes. 
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FIG. 1. Irreversible process of diffusion across a Lorentz slab leading to the equilibration of densities between both reservoirs: 
(a) configuration of the system; (b) initial density profile across the system; (c) time evolution of the densities which decay to 
their equilibrium value in each reservoir. 

The irreversible process of equilibration by the diffusion of particles across the Lorentz slab leads to a global 
production of entropy. Indeed, in the initial situation, the entropy of the whole system is equal to the entropies of 
the ideal gases in both reservoirs if the Lorentz slab is taken to be much smaller than the reservoirs (I <C L) so that 
its entropy is negligible: 



Sir, 



itial 



Sh + Sr. 



with Sl = k B N h In — , and Sr = k B N R In — 
JVL iV R 



with the constant 



7 



ci+ 1 {2irm 1 k B T)i 



according to the Sackur- Tetrode formula [56]. The entropy after equilibration is given by 

2V~f 



Sfinal = k B (N h +N K )]ll 



so that the total entropy increase is 



AS" = S'fina! - ^initial = ^B^L In 



2N L 



2N R 



> k B N L [ 1 



N h + N R 
N L + N R 



k B N R ln 



N^ + N R 



«|1-^I=0 



(84) 
(85) 

(86) 

(87) 
(88) 



2N L J "V 27V R 

because of the inequality ln(l/ir) > 1 — x. Accordingly, the entropy increase (87), which can be rewritten in the form 

(89) 



AS = k B N L In — + k B N R In — > 



is positive in the above process without heat production. 
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We have now to show that the entropy increase (89) is irreversibly produced by diffusion in the Lorentz slab. 
Within the slab, the space and energy distributions f{r, e) of the tracer particles evolve in time according to the 
diffusion equation (49). Since the reservoirs are much larger than the Lorentz slab (L 3> I) the diffusion process is 
quasi-stationary and a linear profile of density f(r, e) is maintained in the slab during the whole relaxation to the 
equilibrium and, this, on each energy shell. The current (51) of particles of energy e is thus transverse to the slab and 
equal to 

J f = -V /r j /l , (90) 

where /r and /l are the densities of particles of energy e in the reservoirs and V = T>(e). These densities relax to 
their equilibrium values according to the equations 

dh _ _dfn _ _Jj_ . , 

dt ~ dt ~ L ■ [ * > 

The entropy of the total system is the sum of the entropies of both reservoirs plus the negligible entropy of the Lorentz 
slab. This total entropy can be expressed in terms of the entropy density s which obeys the balance equation (28) so 
that the time derivative of the total entropy is given by 



dS_ 

dt 



A J dxd t s = A j dx (-d x J s + cr irr ) , (92) 



where A = L d ^ 1 is the area of a section of both reservoirs perpendicular to the direction of the current (see Fig. 
1). The entropy production vanishes in both reservoirs and takes the positive value (79) in the Lorentz slab. Since 
the total system is isolated, the entropy current vanishes at its borders so that the only contribution to the entropy 
variation comes from the irreversible entropy production (79) in the Lorentz slab: 

jc r i-L+l rcc /fl r\2 

— =Adx CT irr = A dx dek B T> 1 xJ> > . (93) 
dt J Jl Jo J 

In a quasi-stationary evolution, the density of particles of energy e is linear in the slab 

/ = h + — j /L (x-L), for L<x<L + l, (94) 



so that the entropy variation (93) becomes 



Tt= k * A L 



de V /r fh In > . (95) 

I JL 



According to Eq. (90), the gradient of densities can be expressed in terms of the current whereupon 

^. = -k B A deJ f ln^. (96) 
dt Jo JL 

Besides, Eq. (91) allows us to write the current as the time derivatives of the densities and we get 

where V = AL is the volume of each reservoir. Eq. (97) can now be integrated over time to obtain 

AS = 5 final - initial = k B V [°° de (/ L In + / R In -j^r-) > , (98) 

Jo V JL + Jr. jl + jr J 

where /l and /r denote the initial distributions given by the Boltzmann distribution (80) with the initial densities 
nh = N^/V and np> = NrJV, while the final distributions are equal to (/l + /r)/2 in both reservoirs. The integral 
over the energy can be performed in Eq. (98) in the same way as the Sackur- Tetrode entropy (81) was obtained. 
Consequently, we recover the entropy increase (89) demonstrating that this entropy is irreversibly produced by diffusion 
in the Lorentz slab. 
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(a) 



heat reservoir 




FIG. 2. Replacement of the Lorentz slab by a movable piston so that work can be extracted from the difference of densities 
if the whole system is in contact with a heat reservoir: (a) initial configuration of the system; (b) final configuration after the 
work W has been extracted. 

Now, the entropy increase (89) corresponds to some work which could otherwise be used if we were to replace the 
Lorentz slab by a movable piston of equal volume, and if the system were put in contact with a heat reservoir (see Fig. 
2). In the initial situation where there is a difference of particle densities between both sides, the piston is submitted 
to a force due to the difference of pressures 

N L k B T N R k B T , N 

-PL.initial = y , and PR,i„itial = — y — • (99) 

This force can produce work if the piston is slowly moved until both pressures equilibrate. We suppose that the piston 
is moved so slowly that the process is isothermal. The equilibrium of pressures is reached when the high-density, 
high-pressure fluid has expanded and compressed the fluid on the other side so that both have the same equilibrium 
density n cq , the temperature T being kept constant. The volumes of both fluids have changed so that 

Nl N r N^ + Nj, N^k B T N R k B T 

n cq = — = -T7- = nT7 > and "L.final = T} = PR,final = , (100) 

with Vl + Vr = 2V, the numbers of particles being here constant in each fluid. The work extracted in this process is 
W = / (P L dV L + P R dV R ) = k B T (n l In ^ + N R In ^f) . (101) 

./initial V V V J 

Since the ratios of volumes are equal to the ratios of densities 

Vl n L V R n R 

— = , and — = , (102) 

V n cq V n cq 

we find that the extracted work is equal to the temperature multiplied by the entropy increase during the previous 
process: 

W = k B T [N L In — + JV R In — ) = T AS . (103) 

Since the initial and final energies of the particles in both fluids are equal by the equation of state E = (d/2)Nk B T 
for ideal gases, the work W has been pumped from the heat reservoir. Accordingly, we conclude that the entropy (89) 



16 



irreversibly produced by diffusion in the Lorentz slab corresponds to some work which is lost during the process of 
diffusion. In this sense, merely by randomizing the velocity directions of the particles, the diffusion process dissipates 
an energy which could otherwise be used, which confirms that diffusion in the Lorentz gas is an irreversible process 
obeying the Second Law of thermodynamics. 



IV. THE IRREVERSIBLE THERMODYNAMICS OF THE MULTIBAKER MODEL OF DIFFUSION 

In this section, we show that the previous considerations for the diffusive Lorentz gas extend to the multibaker 
model of diffusion as well, contrary to the statement by Cohen and Rondoni that "there is no multibaker-analog of 
the concept of local thermodynamic equilibrium" (Ref. [29], p. 125). 



A. The multibaker map as a model of diffusion 



The multibaker model can be considered to be a simplification of the Birkhoff map of the hard-disk Lorentz gas. 
The Birkhoff map of the Lorentz gas is the map induced by the collisional dynamics between the disk scatterers at 
constant energy This map is constructed in terms of two Birkhoff coordinates: (i) the angle < 9 < 2tt of the position 
of collision on the disk, and (ii) the canonically conjugate variable which is the component — 1 < w = smtp < +1 of 
velocity that is tangent to the disk, tp being the angle between the velocity and the normal to the disk at collision. The 
kinetic energy e is conserved during the whole motion so that the trajectories at different energies are similar after a 
rescaling of the time in terms of the speed v(e) = \j2ejra\. Between two collisions, the Birkhoff coordinates £ = (9, zu) 
are constant while the path length A increases linearly with time. Furthermore, the lattice vector t labeling the disk 
on which the collision has occurred is also constant between two collisions. Accordingly, between two collisions, the 
equations of motion become: 



dt 
dX 
~dt 
de 
~dt 
d£ 
~dt 



, (104) 

v(e) = J— , (105) 
V m i 

, (106) 

. (107) 



At collisions, the Birkhoff coordinates, the path length, as well as the lattice vector, £, jump according to the map 

£„ +1 = 0(£„) , (108) 

A„+i = A„+A(£„) , (109) 
= £n , (HO) 

4+1 = In + , (HI) 

where A(£) is the path length of a free flight issued from the impact point £ to the next collision. The lattice 
vector L(£) gives the disk on which the next collision occurs. The time of flight between two collisions is given by 
T(£, e) = A(£)/u(e), which defines the first-return time function or ceiling function. The Birkhoff map together with 
the first-return time function provide an equivalent description of the trajectories of the continuous-time dynamics 
in the Lorentz gas at a given energy. Each phase-space point on a segment of trajectory between two collisions is 
uniquely represented by the coordinates (£, 9, zu, A, e). The hard-disk Lorentz gas has two degrees of freedom so that 
its phase space is four-dimensional. The four original variables of position and momentum (r , p) are here replaced by 
the four continuous variables (6, zu, A, e), the lattice vector £ taking its values in a discrete set. 

In the multibaker map, the Birkhoff map is simplified by considering a baker-type transformation instead of the 
nonlinear transformation <f> induced by the collisional dynamics from disk to disk. Different multibaker models have 
been constructed corresponding to different Lorentz models (see Fig. 3) [10,14-19,23,24,27]. 
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FIG. 3. Schematic representation of different Lorentz gases (I) and their corresponding multibaker model (II) for configu- 
rations of scatterers which are: (a) one-dimensional and periodic; (b) two-dimensional and periodic; (c) one-dimensional and 
disordered; (d) two-dimensional and disordered; (e) with a constant electric field E; (f) with a disordered potential V(x). 
In each example of column (II), the multibaker map acts on the grey squares or rectangles by a baker- type transformation 
combined with translations (depicted by arrows) of the cut pieces to the neighboring squares or rectangles. The multibaker 
(Ha) is defined by Eqs. (112) and (113). 



A dyadic version of the multibaker map is defined by the following mapping which acts on a chain of unit squares 
[16] 




o < x < \ , 

| < x < 1 , 



with the jump vector 



L(x,y) = 



-a 
ha 



< x < \ , 

\<X<1, 



(112) 



(113) 



where a is the spatial distance separating two neighboring squares. The dynamics of the multibaker is time-reversal 
symmetric under the involution I(x,y) — (1 — y, 1 — x) so that <§T X = I o <f>o I. The dyadic multibaker (112) is 
hyperbolic with the positive Lyapunov exponent In 2 and chaotic with the KS entropy /iks = In 2. It is also area- 
preserving as is the Birkhoff map, and the coordinates (0 < x < 1,0 < y < 1) play a similar role as the Birkhoff 
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coordinates (9,w). They also form a pair of canonically conjugate variables. On the other hand, the coordinate 
x corresponds to the unstable direction of the mapping and the coordinate y to its stable direction. This is not 
the case for the Birkhoff coordinates (8, w) which mix the stable and unstable directions of the Birkhoff mapping. 
Therefore, if the linearization of the Birkhoff map of the hard-disk Lorentz gas were possible, a change of variables 
would be required to relate the Birkhoff coordinates (6,w) to the multibaker coordinates (x,y). The Birkhoff map 
of the hard-disk Lorentz gas has singularity lines connected to grazing collisions which are absent in the multibaker 
map. In this regard, the multibaker map constitutes an important simplification of the Birkhoff map of the hard-disk 
dynamics but the analogy is complete as far as all the other properties are concerned, such as the dimensionality of 
the phase space, the diffusivity of the particles moving under the multibaker dynamics, the area-preserving property, 
as well as the positivity of the Lyapunov exponent. It is possible to take into account the energy of the particle in the 
multibaker model by simply keeping this variable as in the Lorentz gas as explained in Refs. [17,18,27]. In analogy 
with the Lorentz gas, energy is conserved. Moreover, we may suppose that the path length between two successive 
intersections in the plane (x, y) of the multibaker map is independent of the coordinates (x, y) and takes the constant 
value A(x,y) = A c . 

In the multibaker map, the phase space is based on a chain of squares. Each square of it is stretched, cut in two 
pieces, and mapped on the neighboring squares. The multibaker is mixing, as is the hard-disk Lorentz gas so that a 
local equilibrium is established in each square of the chain because the dynamics randomizes the coordinates (x,y). 
Moreover, this mapping induces jumps of the particle along the chain of squares. The jumps generate a symmetric 
random walk with diffusion coefficient T> = v(e)a/2, where v(e) is the speed of the particles which may be supposed 
to be the same as in the hard-disk Lorentz gas. 



B. The thermodynamics of the multibaker map 

Except that the dynamics is simplified, the physical interpretation of the multibaker map of diffusion is similar 
to the one of the corresponding Lorentz gas. Thus the thermodynamic considerations of the previous section apply 
mutatis mutandis. Let us consider a gas of Ni particles in a multibaker chain containing N 2 squares. We can impose 
periodic boundary conditions at the ends of the chain for instance. The phase space of the N\ particles is 

M = {(£i,xi,yi,Xi,ei,...,£N 1 ,XNi,yN 1 AN 1 ,eN 1 )\£i + xi < ... < £ Nl + x Nl } , (114) 

with £j e {1, 2, 3, N2}, < Xj < 1,0 < yj < 1,0 < Xj < A c ,0 < ej, with, again, the ordering required by the 
indistinguishability of the particles. In this phase space, the equilibrium microcanonical probability distribution has 
the following density function 

for e < e,- < e + Ae 



F(x 1 ,y 1 ,X 1 ,e 1 ,...,x Nl ,y Nl ,X Nl ,e Nl ) = { V Nl Ae Nl (115) 

otherwise , 

where V is the volume of the domain where the light particles move. This volume can be supposed to depend on the 
physical volume V and on the number N2 of hard disks as in the corresponding Lorentz gas. The equilibrium Gibbs 
entropy (62) calculated for this distribution function is given by 

s - k ^^N^kKx' (116) 

where Ax Ay AX Ae is the volume of the cells partitioning the phase space. The Gibbs relation (65) is satisfied with a 
vanishing inverse temperature, the chemical potential of the light particles given by 

Mi NjAxAyAX 

— = k B ln , (117) 



and another one for the scatterers 



H2 N x k B ( dv 



T V \ 8N< 



(118) 



Hence, the considerations of the previous section follow, showing that the multibaker map is also consistent with the 
irreversible thermodynamics of diffusion processes, contrary to the claims by Cohen and Rondoni [29]. 
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In Refs. [17,18,20,23-27,53,54], the entropy production associated with the transport processes sustained by multi- 
baker models was derived and the expressions expected from irreversible thermodynamics have been obtained. More 
recently, the ab initio derivation of entropy production of diffusion has been extended to the continuous-time dynamics 
of the Lorentz gas as well as to many-body systems composed of a tracer particle moving in a fluid of other particles 
[55]. These works show that the thermodynamically expected entropy production can be successfully derived from 
dynamics in the Lorentz gases and multibaker maps. 



V. CONCLUSIONS 



This paper provides counter-arguments refuting the doubts and criticisms expressed by Cohen and Rondoni con- 
cerning the thermodynamic considerations of Lorentz gases and multibaker maps. We have shown that their claims 
cannot be taken uncritically, and do not stand up under a rather direct consideration of the irreversible thermody- 
namics of isothermal mixing of distinct chemical species. We hope that this exchange of ideas will be helpful for 
clarifying the thermodynamics of the model systems we study. We feel that, at best, the arguments of Cohen and 
Rondoni should be considered as restatements of the facts that a study of simple systems may not reveal all of the 
features and difficulties that one encounters in a study of complicated ones. Yet, whatever the specificities of simple 
systems may be, we have seen that they remain within the realm of irreversible thermodynamics. The present work 
thus shows that irreversible thermodynamics indeed has a universal applicability in the class of spatially-extended 
many-particle systems. 

We should point out that the goal of the line of research presented in the papers criticized by Cohen and Rondoni is 
not the treatment of entropy production in Lorentz gases or multi-baker maps. Rather, these papers represent steps 
on the way to a treatment of entropy production in classical fluids. As is often the case in statistical physics, the 
study of simple models allows workers to develop insights into the ideas and techniques needed to treat more realistic 
systems. From this point of view the criticisms of Cohen and Rondoni, while mistaken, as we argue here, would also 
be irrelevant if our line of research shows that similar methods can be useful in the study of fluid systems where all of 
the particles can interact with each other. In fact, we have shown in a previous paper [55], that the methods under 
discussion apply equally well to a periodic, fluid system with a tracer particle. In this case, we considered a fully 
interacting, periodic N + 1 particle system, in a microcanonical ensemble. For this system the presence of interacting, 
moving particles does not change the arguments in any important way. 

Cohen and Rondoni have also raised more technical objections to our use of Sinai- Ruellc-Bowen measures to calculate 
the rate of entropy production as described in Refs. [14,20,53-55] which also require a response. This will be supplied 
in a future publication. 
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